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ON  THE  DETECTION  OF  ULTRA  WIDEBAND  RADAR  SIGNALS 


1.  INTRODUCTION 

An  ultrawideband  (UWB)  radar  is  a  radar  that  has  an  instantaneous  bandwidth  which  is  relatively 
large.  An  impulse  radar,  which  is  one  type  of  UWB  radar,  is  a  radar  that  has  a  time-bandwidth  product 
nearly  equal  to  1.  An  introduction  to  impulse  radar  can  be  found  in  Ref.  1. 

The  primary  use  of  UWB  radar  is  in  applications  in  which  conventional  radars  do  not  perform 
satisfactorily.  Fowler  et  al.  [2]  concluded,  however,  that  impulse  radars 

•  show  no  evidence  of  phenomenological  advantages, 

•  have  no  unique  advantages  in  terms  of  anti-stealth,  and 

•  possess  no  special  low  probability  of  intercept  (LPI)  advantages. 

But  Fowler  et  al.  also  concluded  that  UWB  radars  may  prove  to  be  more  effective  over  conventional 
radars  in  certain  short-range  applications  which  require  high  range  resolution  (HRR)  such  as  terrain 
profiling.  Reference  3  lists  target  identification,  detection  of  specific  air-turbulence  patterns,  and 
detection  of  low  radar  cross  section  targets  as  other  potential  HRR  applications. 

An  obstacle  to  implementing  a  true  wideband  radar  is  the  tremendous  amount  of  data  which  needs 
to  be  processed  in  real  time.  Another  problem  is  in  generating  the  energy,  within  a  period  typically  on 
the  order  of  picoseconds,  necessary  to  detect  or  discriminate  a  target. 

This  report  examines  general  signal  processing  aspects  of  UWB  radar  from  a  detection  standpoint. 
In  Section  2,  we  model  the  process  and  derive  a  detector  structure  for  a  point  target  in  white  Gaussian 
noise.  Section  3  derives  a  new  generalized  moving  target  indicator  (MTI)  detector  design  for  cases  of 
correlated  Gaussian  noise  or  clutter,  and  extends  the  detector  to  a  generalized  Doppler  filter  bank. 
Additionally,  the  associated  effects  of  bandwidth  on  performance  in  the  presence  of  clutter  are  examined. 
Section  4  considers  range  spread  targets  and  derives  three  fundamental  detector  structures.  Section  5 
discusses  future  work  areas  and  summarizes  the  report. 

It  should  be  mentioned  that  much  of  the  work  is  based  on  extending,  when  necessary,  conventional 
radar  processing  theory  taken  from  Van  Trees  [4,5}. 

2.  POINT  TARGET  DETECTION 

We  will  now  derive  the  optimal  detection  statistic  for  an  UWB  signal  assuming  a  point  target.  We 
assume  the  (single  pulse)  transmitted  waveform  is  given  by 


exp(/oy)]’ 
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where  Et  is  a  constant  representing  the  transmitted  energy  of  the  signal  and  f(t)  is  the  complex  envelope 
of  the  transmitted  signal  wc,  the  carrier  frequency.  It  is  assumed  that 


}  \m\2dt  =  i. 

We  assume  that  the  propagation  and  reflection  processes  are  linear  and  frequency-independent,  and  that 
the  target  is  at  a  distance  Ra  at  t  —  0  and  is  moving  with  constant  radial  velocity  v.  Under  these 
circumstances,  it  is  straightforward  to  show  [4,  p.  241]  that  the  received  signal  is 


\f2Re\jE^Bf(t  -  r(t))  exp (jac(t  -  r(/)))]. 


where 


r(r)  = 


2  RJc 
1  -  vie 


2vtlc 
1  -  vie 


is  the  round-trip  time  delay,  b  represents  an  unknown  complex  amplitude,  which  is  often  assumed  to  have 
a  zero  mean  Gaussian  distribution.  We  can  reduce  exp(/wf(f  -  r(r)))  by  absorbing  the  first  term  of 

2vw  Ic 

r(r)  into  b  and  defining  w .  =  -  as 

1  ~  v/c 

exp(/'wc(r  -  r(r))  «exp (j(uc  -*■  wd)t). 

To  simplify  f{t  -  r(t)),  the  approximation  r(r)  =  2 RJc  -  Ivtlc  is  made  for  conventional  radars. 
However,  the  actual  time  the  signal  return  begins  (which  we  will  refer  to  as  tb )  assuming /begins  at  t 
=  0,  is  when  r(fb)  =  tb,  which  simplifies  to 


2RJe_ 

1  +  v/c ' 


whereas  the  latter  approximation  yields 


,  _  2  RJc 

b  1  +  2  vie 

Now  suppose  we  are  trying  to  detect  a  high-speed  target  with  Rp  =  20,000  m  and  u  =  600  m/s.  The 
approximation  yields  rb  ~  1.333328E  -  4  whereas  the  actual  time  is  rb  =  1  333301E  -  4.  a  difference 
of  2.7-10.  Such  a  difference  may  seem  insignificant  until  it  is  compared  with  the  extent  of  a  typical 
pulse.  For  example,  assuming  a  bandwidth  of  5  GHz,  the  pulse  width  is  on  the  order  of  a  few  hundred 
ps,  or  of  the  same  order  as  the  time  difference  between  the  approximated  value  and  the  actual  value. 
Thus,  we  could  be  in  range  error  by  as  much  as  one  range  cell  width  when  using  this  approximation.  We 
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can  now  accurately  approximate  fit  -  r(t))  by  fit  -  r)  where  r  =  - - - Note  that  we  have 

1  +  vie 

removed  the  dependence  of  r(t)  on  t  by  replacing  r(r)  by  the  epoch  tb  of  the  received  pulse.  This  is 
valid  as  long  as  the  compression  or  stretching  of  the  time  scale  caused  by  the  dependence  on  time  is  not 
significant.  It  is  shown  in  Ref.  4  (p.  241)  that  this  is  the  case  whenever  the  time  bandwidth  product  of 
the  received  signal  is  <  c/(2i>);  this  will  be  assumed  for  the  UWB  signals  considered  here. 

We  assume  that  one  of  two  hypotheses  is  received.  The  hypotheses  are  signals  present  with  additive 
noise  n(t)  under  Hj  and  only  additive  noise  under  H0.  Under  hypothesis  H,  we  receive 


r(t)  =  JirJJeJ~ B fit  -  r)  exp(/7(wf  +  ud))  +  n(t)  exp(/wcf)]. 


(1) 


while  under  H0 


r(r)  =  i/2 Re[nit)  exp(/«cr)]  (2) 

is  received.  ud  is  an  added  Doppler  frequency  shift.  We  now  make  the  additional  assumption  that  the 
signal  is  band-limited.  Under  these  circumstances,  there  exist  the  transformations 

rcit)  =  i'Jlrit)  cosuf)^ 
rsiO  =  i'Jlrit)  sina)^)^, 
rit )  =  J2rcit )  cosw^  +  Jlrft)  sinw/ 

where  i  )u,  represents  a  low-pass  operation.  Note  these  operations  are  invertible,  hence  by  the 
factorization  theorem  [5]  it  can  be  shown  that  the  problem  is  equivalent  to  discriminating  between 


rit)  =  ^E^bfit  -  t)  exp(jwdt)  +  nit)  under  H,  and 

rit)  =  nit)  under  H0.  (3) 

This  problem  has  the  solution  [4,  p.  245] 


", 

I  1 1'  exP i-jujOdt  | 2  <  y  , 

"0 

where  Tp  is  the  pulse  repetition  interval  and  7  is  a  threshold  associated  with  the  probability  of  false  alarm. 

The  test  is  uniformly  most  powerful  for  the  case  where  b  =  $  exp (J4>),  f  is  unknown  and  0  is  a 
uniformly  distributed  random  variable  on  [0,2tt].  To  prove  this,  examine  Ref.  5  (p.  339).  The  test  is 

optimal  in  the  case  where  5  is  assumed  Gaussian  distributed  with  zero  mean,  which  is  the  case  when  f 
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is  Rayleigh  distributed  and  <f>  is  uniformly  distributed  on  [0,  2x][7],  For  the  case  of  n  pulses,  this 
becomes 


|  If  '■(O  E  /*('  ■  Ti)  exp(-;aj/)dr  |2  <  y 

«o 


or 


I  E  (n  ^  f  (')/*('  “  r()  txp(-judt)dt  f  y, 

,  -  o  * 0 

"o 


where  r  =  ivTp),c  +  ,r  j  =  0, n  -  1.  This  detector  assumes  knowledge  of 

1  l+u/c 

_  2vuclc  whjch^  in  practice>  js  typically  not  known.  An  implementable  detector  structure  is 
1+u/c 

developed  in  Section  3  that  does  not  assume  knowledge  of  ud  and  is  extended  from  a  white  noise 
assumption  to  correlated  noise  or  clutter. 

3.  DETECTION  IN  CLUTTER 

In  the  derivation  of  Eq.  (4),  the  target  is  not  assumed  to  lie  in  the  same  range  cell  from  pulse  to 
pulse.  The  associated  delays  r,  are  dependent  upon  the  target  velocity  in  contrast  to  a  conventional  radar 
where  r,  =  iTp  is  nearly  velocity  independent.  In  the  presence  of  clutter  or  correlated  noise,  this 
difference  has  important  consequences.  In  the  conventional  case,  the  clutter  within  the  target  range  cell 
can  be  cancelled  using  only  data  from  returns  of  the  target  at  the  same  range  cell.  Hence,  the  same  data 
can  be  used  to  both  cancel  clutter  and  integrate  the  target  energy  from  successive  pulses.  In  the  UWB 
case,  however,  the  target  may  move  through  range  cells  between  successive  pulses  so  that  the  range  cells 
that  contain  the  target  may  not  be  sufficient  to  cancel  the  clutter.  That  is,  the  degree  of  clutter  correlation 
between  successive  target  returns  that  lie  in  different  range  cells  may  be  much  less  than  in  the  case  of 
conventional  radars.  In  general,  we  would  require  additional  data  that  do  not  contain  the  target  to  fully 
cancel  the  clutter  and  obtain  reasonable  detection  performance. 
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We  refer  again  to  the  hypothesis  testing  problem  of  Eqs.  (1)  and  (2)  where  we  model  the  cluttern(r) 
by  a  correlated  zero-mean  nonwhite  Gaussian  noise  process.  The  optimal  detector  is  found  by  using 
standard  techniques 


Jo"7'  ?(t)g  *(0  exp(-judt)dt 


", 


> 


(5) 


where  g(t)  is  the  pseudosignal  formed  by  solving 


E  /('  ■  Ti)  =  [  "Tp  Kn(t,u)g{u)du 

i  »  0  J0 

and  Kn(t,u)  =  E(n{t)n  *(«)).  The  pseudosignal  would  take  into  account  correlation  in  range  as  well 
as  time  correlation  from  pulse  to  pul  e. 

Generally,  the  Doppler  ud  is  unknown.  A  bank  of  filters  of  the  form  of  Eqs.  (4)  or  (5)  may  be  used 
such  that  each  filter  is  matched  to  a  different  Doppler  shift.  When  the  maximum  of  these  filters  is 
compared  with  a  threshold,  the  resulting  test  is  often  referred  to  as  a  generalized  likelihood  ratio  test 
(although  an  infinite  number  of  filters  is  technically  required).  In  this  case,  we  would  require  a  number 
of  integrators  to  perform  the  computation  of  either  Eqs.  (4)  or  (5). 

Let  us  denote  the  output  of  a  filter  matched  to  /(r)  by 

u(t)  =  jr(r)f(T  -  t)dr.  (6) 

Suppose  that  nearly  all  the  energy  of /(/)  is  concentrated  within  the  time  interval  [0,  7^-1 .  If  wdT  varies 
relatively  little  during  this  interval,  we  can  make  the  approximation  [8] 

u(t)  *  exp(judt)^r(T)f(T  -  t)  exp(-jwdT)dT. 

For  the  case  of  white  Gaussian  noise,  the  detector  structure  of  Eq.  (4)  then  reduces  to  the  optimal  solution 


", 

n  - 1  j  > 

|  Y,  exp(-ywd7(.)«(r()  |  <  y. 

i-O 

"o 


Although  Eq.  (6)  is  a  sufficient  statistic  for  the  detection  problem  in  white  Gaussian  noise,  it  is  not  in 
general  sufficient  for  arbitrary  correlated  noise.  However,  in  the  case  of  a  conventional  radar 
(7,  =  iTp),  a  complex  weighting  a,  is  typically  substituted  for  the  values  expf-y'c^r,)  in  Eq.  (7).  The 
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use  of  an  appropriate  weighting  is  well  known  to  eliminate  the  effects  of  most  clutter.  The  resulting 
detector  is  then 

", 

Ip  ,  ,  |2  >  (8) 

E  <  y- 

i  =  0 

"o 


This  detector  for  appropriate  a,  is  referred  to  as  an  MTI.  A  bank  of  such  filters  is  a  Doppler  filter  bank. 

In  the  case  of  UWB  operating  in  clutter,  the  detector  structure  of  Eq.  (8)  may  not  be  sufficient  to 
cancel  the  clutter.  This  is  due  to  the  loss  of  correlation  between  the  m(t,)  which  may  occur  if  the  target 
moves  between  range  cells  on  successive  pulse  repetition  intervals.  A  generalization  of  Eq.  (8)  is  thus 


", 

E  E  aiju(Ti +  tp(j  ~  0)  |  <  7 

(=0  j  =  0 

"o 


(9) 


where  n  is  the  number  of  pulses.  We  define  a  detector  of  the  form  described  in  Eq.  (9)  to  be  a 
generalized  MTI  (GMTI),  and  the  n  by  n  matrix  A  of  elements  atJ  to  be  the  GMTI  coefficients.  We  can 

(y  .) 

extend  the  GMTI  to  a  generalized  Doppler  filter  bank  as  follows.  Let  atJ  be  a  set  of  Doppler 
coefficients  matched  to  a  target  with  associated  Doppler  We  can  compare  each  filter  with  a  threshold 
or  declare  a  detection  whenever 


Performance  in  the  white  Gaussian  noise  case  is  well  known.  We  can  relate  the  probability  of 
detection  Pd  with  the  probability  of  false  alarm  Pf\ 4,  p.  246J 

pf=py*  (10) 

where  A  =  E/N0,  E  is  the  expected  value  of  the  received  signal  energy  and  N0  is  the  (complex)  noise 

spectral  density.  In  the  case  of  correlated  Gaussian  noise,  Eq.  (10)  holds  although  A  =  E  |/(f)g  *  {t)dt. 
In  either  case,  performance  is  limited  by  the  spectral  density  of  the  thermal  noise  and  the  received  signal 
energy.  This  is  also  the  case  for  conventional  radar. 
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3.1  Computation  of  GMTI  Coefficients 

In  this  section  we  show  how  to  compute  the  GMTI  coefficients  of  Eq.  (9).  Consider  the  vector  x 
where  xa  =  u(t0),  =  u(r0  +  7»,  =  u(T0  +  (n- l)Tp) ,  xn  =  uO,-7»,  xn =  u(t0), 

m(t0),  x^_x  =  u(t I  +  (n  -  2)TP)  continuing  through  xnM.  In  general, xjn+;  =  u(t,  +  Tp(j  -  i)) 
where  0  £  i  £  n  -  1  and  0  <y  <  n  -  1  (see  Fig.  1).  We  see  that  the  entries 
will  be  correlated  because  of  the  temporal  correlation  of  clutter,  while  Xj,xn^, ...  will  be 

correlated  because  of  the  spatial  correlation  of  clutter.  Other  cases  represent  combinations  of  both  time 
and  spatial  correlation.  Assuming  stationarity,  we  can  represent  the  covariance  matrix  R  of  x  by 


R  =  E(xx11)  = 


R,  Rh  Rs  R  1 

‘  s>  s»-2  s» 

R .  R.  R.  R . 


7?  «  7?  «  •••  7?.  7? 

'»-2  ■'»-!  '  Jl 

7?  h  7?  h  •  7?  7? 

s«-l  s»-2  ' 


where  /?,  represents  the  covariance  matrix  associated  with  any  row  of  Fig.  1  and  Rs  is  the  cross 
covariance  matrix  associated  with  any  two  rows  of  Fig.  1  that  are  separated  by  i  rows.  xH  denotes  the 
Hermitian  transpose  of  x.  Hence,  R,  is  a  clutter  covariance  matrix  that  may  have  relatively  large  values, 
while  the  matrices  Rs  will  typically  diminish  with  increasing  /  as  the  clutter  decorrelates  spatially. 


Fig.  1  —  Matched  filter  output  array  for  a  target  with  range  decreasing  in  time. 
Target  occupancy  is  XQ,  ...,  Xn7^. 


Assuming  Gaussian  noise  n(t),  the  vector  x  is  easily  shown  to  be  Gaussian.  We  can  formulate  the 
new  detection  problem  as 


y  =  x  +  s  exp (j<t>)  under  H,  and 
y  =  x  under  H0, 
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where  s  is  the  expected  value  of  x  under  H,  and  <f>  is  a  random  uniformly  distributed  phase.  It  is  known 
[9]  that  the  optimal  solution  is 


", 

n2- I 

> 

E  v,| 

< 

i  =  0 

"0 

(11) 

' 


where  c  =  R  'V.  To  illustrate  this,  we  assume  for  simplicity  that  n  =  2,  s  =  (1,0,0, 1)H  and 


1  p  0  0 
p  1  0  0 

R  = 

0  0  1  p 
0  0  p  1 

The  optimal  coefficient  vector  is  then  given  by  (l,-p,-p,l).  As  p  approaches  one,  the  optimal  solution 
becomes 

", 

> 

\Xo  ~  X\  “  *2  +  *3  i  <  y 

"o 


which  essentially  removes  the  clutter  and  integrates  the  signal.  Note  that  the  signal  only  lies  in  two 
components.  A  similar  result  could  easily  be  derived  to  allow  implementation  with  a  transversal  filter. 
Hence,  in  cases  such  as  this,  Eq.  (8)  would  not  be  adequate.  This  situation  is  extreme  in  that  the  target 
returns  are  separated  by  at  least  a  range  cell  between  pulse  repetition  intervals.  In  cases  where  this  does 
not  occur,  the  complex  representation  of  Eq.  (3)  is  very  useful,  as  the  optimal  detector  of  Eq.  (11)  will 
incorporate  the  phase  change  of  the  signal  relative  to  the  carrier  frequency,  allowing  relatively  slow 
targets  to  be  detected. 

3.2  Spatial  Correlation 

Increasing  bandwidth  or  resolution  may  have  the  effect  of  spatially  correlating  range  cells  compared 
with  a  narrowband  radar.  The  effect  that  spatial  correlation  might  have  on  the  signal-to-noise  ratio  can 
be  determined  via  measurements  made  to  estimate  the  covariance  matrix  R.  Let  R ,  be  the  covariance 
matrix  where  entries  due  to  spatial  correlation  are  zero.  It  is  known  [9J  that  the  optimal  improvement  in 
signal-to-noise  ratio  is  given  by 

—  |sTR "'s  *1  ,  (12) 

7 

where  s  is  any  signal  vector  corresponding  to  Eq.  (1 1)  and  7  is  the  received  signal-to-noise  ratio.  We 
simply  mention  here  that  a  measure  of  relative  performance  is  given  by 

x  - 

|s tR,’is  *  | 
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3.3  Clutter  Resolution 

Increasing  the  resolution  will  have  an  impact  on  a  radar’s  ability  to  resolve  a  limited  number  of 
clutter  scatterers.  The  problem  will  be  examined  by  considering  the  processing  of  a  single  transmitted 
pulse  without  Doppler  processing.  The  environment  is  assumed  to  consist  of  a  target  between  two  clutter 
scatterers.  The  detection  performance  is  investigated  as  a  function  of  both  the  bandwidth  and  the  distance 
a  target  is  between  the  two  clutter  scatterers.  The  general  resolution  problem  was  studied  in  detail  by 
Van  Trees  [4,  pp.  323-335],  although  we  are  interested  in  the  effects  of  bandwidth  on  clutter  resolution. 

We  fix  the  distance  between  the  two  clutter  scatterers  at  1  m  for  the  examples  to  follow  simply  to 
illustrate  the  effects  of  bandwidth  on  clutter  resolution  by  conventional  and  optimal  processing.  We 
assume  that  a  point  target  is  located  at  a  distance  d  from  cluttei  scatterer  number  1  as  indicated  in  Fig. 
2.  Let  f(t)  be  the  received  signal  and  for  simplicity  zero  target  Doppler  is  assumed.  The  received 
signals  from  scatterers  1  and  2  are,  respectively,  bxf(t  -  7,)  expOlO  and  b^fit  -  r2)  exp(/co2r), 
where  t,  =  -Idle,  and  r2  =  2(1  -  d)lc.  Hence  we  formulate  the  hypothesis 

2 

r(t)  =  +  y /i^" Y,  "  Ti )  exp 0,0  +  n(t)  under  Hj  and 

i  =  1 

2 

r{t)  =  £  b/(t  -  r,)  expO,0  ■f  n(t)  under  H0. 

1=0 


T  •  Sr 


cc 

LU 

H- 

LU 

5 


TARGET 


Fig.  2  —  Depiction  of  clutter  scenario.  S ,  and  S2 
arc  clutter  scatterers  separated  by  I  meter. 
Target  is  at  d  meters  from  Sv 


The  bl  are  assumed  to  be  zero  mean  complex  Gaussian  random  variables  such  that E(bjb*)  =  2  a^5f. 
where  5tJ  =  1  for  i  =  j  and  0  otherwise.  n(t )  is  zero  mean  white  Gaussian  noise  with  covariance 

function  E(n(t)n  *(«))  =  N05(t  -  u) .  We  first  assume  the  detector  is  of  the  conventional  type,  and 
later  assume  the  detector  is  optimal.  In  either  case,  the  performance  of  this  detector  can  be  derived  as 
|4,  p.  251] 


r/u 


=  p 
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where  the  performance  exponent  A  =  Ac  for  the  conventional  receiver  and  A  =  A, 
receiver.  The  conventional  detector  is  given  by 

"i 

|  Jr(0/*(/)df  |  <  7 


and  the  performance  coefficient  [4,  p.  327]  is  given  by 


A, 


Eo 

No 


1  +  ^  +  ~ 0(t2,o)2 ) 

"0  "o 


where  £.  =  2 Eto* ,  and  6 (r(, )  is  the  ambiguity  function  of  /(r), 

0(t,w)  =  |  |/(f  “  ^)/*(*  +  ~)  exp (jut)dt  f. 

The  conventional  receiver  is  not  optimal  due  to  the  inclusion  of  the  clutter  scatterers. 
[4,  p.  330]  that  the  optimal  receiver  is 


Hi 

|  jr(r)«  |  <  7 

"o 

where 

2 

g( 0  =  8  JO)  +  Y,  SjO  -  Ti)  expO'w.r) 

i  =  1 

and  the  coefficients  g,  are  computed  as  follows.  Define: 


// 


A 

P 


Pd 


f(t  -  Tj)  exp(yw,0 
f{t  -  t2)  exp(yw2/) 

X  0 
0  T2 

J  f'iOf'itfdt 


for  the  optimal 


(13) 


It  can  be  shown 
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Then  gQ  =  —  and 


*  ■  U' I  ■ -i?  [' • 


Additionally,  the  optimal  exponent  Aa  can  be  derived  as 


7T  r  -i 

.  1  1  H  ,  1  .  .  . 

’  7T  1  ‘  W*'  '  ^ 

o  L  o  ’o 


_  A  A 

We  assume  w,  =  0  and  E  =  CIN  for  i  =  1,2,  and  define  R.  =  - —  and  Rn  =  — — .  Rc  and  Ra 

EJN0  EJN0 

represent  normalized  signal-to-noise  ratios  after  conventional  and  optimal  filtering,  respectively.  In  this 
case,  Eq.  (13)  can  be  rewritten  as 


1  +  —0(7-.,  0)  +  —  0(7,,  0) 

N0  ‘  K 


and  Eq.  (14)  as 


*o  =  1  + 


1  +  — —  -0  7-.-  r,,0 

CIN  1  2 


f  (e(TV°)  +  6(T2’°) 


-2  0(r,  -  r2, 0)0(7,, O)0(r2,o|  , 

> 


where  7,  and  t2  are  parameterized  by  d  as  in  Fig.  2.  Hence,  the  performance  of  both  the  conventional 
and  optimal  receivers  depends  on  the  ambiguity  function  of  /(f).  For  example,  let  /(f)  be  a  Gaussian 
pulse 


/(O  =  [-4  1  exP  ~~~2 
7 tT2  2  T2 


It  is  known  that  the  ambiguity  function  is  |4,  p.  283] 


0(7, w)  =  exp  -1  1-  +  TWj  j , 
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where  7  is  a  constant  proportional  to  the  useful  duration  of  the  pulse.  It  is  easily  shown  that  the  power 

2  2 

spectrum  of  f(t)  is  proportional  to  e'T'a"  and,  hence,  the  3  dB  bandwidth  uJdB  =  radians  per 

second. 

Assuming  a  Gaussian  pulse,  we  plot  in  Figs.  3{a)-3(d)  Rc  vs  d,  and  in  Figs.  4(a)-4(d),  Ra  vs  d,  for 
various  C/N  ratios  and  u3dB  bandwidths.  It  is  seen  that  at  100  MHz,  detection  performance  is  limited 
by  the  clutter  residue,  while  at  higher  frequencies,  much  of  the  region  is  free  of  the  clutter  from  a 
detection  standpoint.  Also  note  in  cases  such  as  o>3dB  =  1  GHz  the  optimal  solution  substantially 
improves  subclutter  visibility  of  the  target  echo  over  the  conventional  processor. 


0  0)  0  2  03  04  05  0  G  07  08  09  1 

d  (m) 

Fig.  3(a)  —  Rc  vs  d,  frequency  =  100  MHz, 
conventional  processor 


Fig.  3(b)  —  Rc  vs  d,  frequency  =  1  GHz, 
conventional  processor 


Fig.  3(c)  —  Rr  vs  d,  frequency  =  5  GHz,  Fig.  3(d)  —  Rc  vs  d,  frequency  =  10  GHz, 

conventional  processor  conventional  processor 
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Fig.  4(a)  —  Ra  vs  d,  frequency  = 
optimal  processor 

100  MHz, 

j _ i  i _ I _ |_ 


0  0.1  0.2  0.3  04  0.5  0.6  07  C.8  09  1 

d(m) 


Fig.  4(b)  —  Ra  vs  d,  frequency  =  1  GHz, 
optimal  processor 


Fig.  4(c)  —  Ra  vs  d,  frequency  =  5  GHz, 
optimal  processor 


Fig.  4(d)  —  Ra  vs  d,  frequency  =  10  GHz, 
optimal  processor 


Implementation  of  the  optimal  processor  requires  knowledge  of  the  positions  of  the  clutter  scatterers 
and  their  relative  strengths.  Such  processing,  however,  may  be  achieved  adaptively  by  estimating  the 
clutter  strengths  and  positions  assuming  a  weak  target  signal.  Such  a  processor  may  be  useful  to  reduce 
the  effects  of  clutter  when  postprocessing  is  not  sufficient  to  fully  cancel  clutter. 

4.  RANGE  SPREAD  TARGETS 

Due  to  the  fine  resolution  of  HRR  radars,  target  returns  will  typically  occupy  more  than  one  range 
cell.  We  consider  three  cases.  In  the  first  case,  the  target  amplitude  profile  is  assumed  known  up  to  a 
proportionality  constant.  An  optimal  detector  is  derived  which  is  an  extension  of  Eq.  (7).  In  the 
remaining  two  cases,  the  target  profile  is  not  known  exactly,  although  certain  second  order  statistics  are 
assumed  known  regarding  the  target.  For  the  second  case,  we  assume  that  the  target  scatterers 
individually  resolve  and  model  the  return  as  a  sum  of  resolved  point  scatterers  with  statistically 
independent  Gaussian  amplitudes.  For  the  final  case,  the  target  is  modeled  by  a  weighted  integral  of  the 
transmitted  signal  where  the  weighting  function  is  assumed  zero  mean  Gaussian  with  known  covariance 
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function.  The  returns  from  different  ranges  are  assumed  statistically  independent.  An  optimal  detector 
is  found  via  a  low  energy  coherence  (LEC)  assumption  [4,  p.  431].  The  LEC  case  requires  that  the 
eigenvalues  X;  in  the  Karhunen-Lo&ve  expansion  of  the  signal  s(t )  are  small  compared  to  the  spectral 
density  of  the  white  noise  level.  Since  the  signal  energy  is  spread  across  a  large  bandwidth  in  UWB 
radar,  this  is  a  reasonable  assumption. 

For  simplicity,  we  will  assume  the  target  is  embedded  in  white  Gaussian  noise  and  consider  the 
hypothesis  test  in  Eq.  (3).  Suppose  f(t )  can  be  written  as 

/( 0  =  j*(X)s(f  -  \)d\, 

where  b(\)  is  a  known  complex  function  proportional  to  the  amplitudes  of  the  target  scatterers  and.v(/) 
is  the  transmitted  signal.  Let  us  denote  the  output  of  a  filter  matched  to  s(t)  by 

us(t)  =  -  t)dT.  07) 

The  optimal  receiver  of  Eq.  (7)  is  given  by 


Note  that  if  the  individual  target  scatterers  are  resolved,  the  term  f  b  *  (\)us{Tj  +  \)d\  essentially  is 
a  coherent  addition  of  the  individual  target  scatterers  energies.  If  the  scatterers  do  not  resolve,  the 
amplitudes  will  typically  be  such  that  some  cancellation  occurs  between  the  scatterers  resulting  in  a  loss 
of  energy.  This  effect  was  noted  by  Farina  |3]  who  showed  that  for  a  moderate  number  of  scatterers, 
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significant  improvements  in  the  detection  performance  are  possible  by  using  an  UWB  over  a  low- 
resolution  radar. 

In  the  second  case,  we  will  assume  the  target  scatterers  are  resolved  and  that  their  amplitudes  are 
independent  Gaussian  random  variables.  That  is. 


m  =  E  M('  -  \) 

i  -  0 


where  the  bi  are  zero  mean  Gaussian  with  E(btb*)  =  2a]  b^.  If  the  X,  are  separated  enough  so  that 
s(t  -  X()  can  be  assumed  orthonormal,  the  signal  covariance  matrix  is  separable,  and  the  optimal 
solution  is  [6,  p.  352] 


T->!  2a?  |  ,  |  > 

E  — -  <  7 

.  ■  o  2a-  +  N0  J 

" o 


For  n  multiple  pulses,  E"=o  “  T,)>  we  will  assume  that  the  scatterers  maintain  the  same  geometric 
phase  relation  from  pulse  to  pulse  resulting  in  the  detection  statistic 


K  -  1 


E 


2a' 


2a'  +  Na 


exp (-judTf)  |f(r)5(r  -  X, 


(19) 


Again  we  see  that  the  individual  scatterer  energy  is  added  to  form  the  detection  statistic. 

Generally  the  function  b(\)  is  highly  target  dependent,  or  may  not  be  known.  Suppose  b(\)  is 
modeled  as  a  sample  function  from  a  zero  mean  complex  Gaussian  process  5(k)  such  that 

f{t)  =  p(X)s(r  -  \)d\ 


and  furthermore  assume  the  covariance  function  of  b(\)  is  of  the  form 

^(X,^)  =  5 (Xj  -  XjWX,) 

where  5(X)  =  E(\f>(X)b  * (X) | 2 ,  and  5(*)  represents  the  impulse  function.  The  optimal  solution  in 
general  is  difficult  to  implement  because  there  can  be  an  infinite  number  of  eigenvalues/eigenfunctions 
in  the  Karhunen-Lo^ve  expansion  of  ^(XpXj).  However  if  we  make  the  LEC  assumption  that  the 
eigenvalues  in  a  Karhunen-Lo£ve  expansion  are  much  less  than  the  spectral  density  of  the  white  Gaussian 
noise,  approximations  can  be  used  which  simplify  the  detector  structure.  For  the  UWB  case,  this  is 
typically  the  case  for  weak  targets  due  the  wide  bandwidth  of  the  signal.  For  strong  targets,  the  detector 
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structure,  although  suboptimum,  may  still  perform  sufficiently  well,  since  the  strength  of  the  target  may 
suffice  to  overcome  the  deficiencies  of  the  suboptimum  detector.  For  a  single  pulse,  the  optimal  detector 
is  [4,  p.  431 1 


JS(\)K(X)|2d\  <  7, 

«o 

which  is  easily  extended  to  n  multiple  pulses  £  "=  o  f(t  -  r(). 


n  - 1 

S{\)  |  Y,  exp (-j<j)dTi)us(Ti  +  X)  |  dk 

i  =  0 


We  have  compiled  three  optimal  detector  structures  for  range  spread  targets.  The  actual  detector 
structure  which  is  most  applicable  to  UWB  depends  on  the  particular  targets  of  interest  to  the  radar 
designer  and  feasibility  of  implementation. 

5.  CONCLUSIONS 

The  purpose  of  this  report  is  to  examine  general  aspects  of  the  detection  of  UWB  signals  and  cite 
areas  for  further  study.  We  determined  detectors  assuming  white  Gaussian  noise  and  correlated  noise, 
found  a  detection  structure  which  generalizes  a  conventional  Doppler  filter  bank,  and  illustrated  a  method 
to  compute  the  GMTI  coefficients.  We  also  examined  the  effects  of  increasing  bandwidth  on  clutter 
resolution  and  compared  optimal  and  suboptimal  detectors  for  single  pulse  echos.  We  found  that  the 
optimal  processor  can  result  in  substantial  improvement  in  subclutter  visibility  of  the  target  echo.  We  also 
examined  range  spread  or  distributed  targets.  Cases  of  completely  known  and  incomplete  statistical 
knowledge  of  the  target  profile  were  treated. 

The  determination  of  the  parameters  of  many  of  the  detection  structures  discussed,  such  as  optimal 
weights  of  a  filter  bank,  has  been  shown  to  depend  on  the  statistical  characterization  of  the  signal  return 
which,  in  turn,  depends  on  many  factors.  Signal  characteristics  (such  as  form  and  bandwidth),  target 
characteristics  (such  as  form  and  velocity),  and  the  environment  (which,  along  with  various  radar 
parameters,  determines  the  degree  of  clutter  and  spatial  correlation  and  the  signal  propagation  paths)  are 
all  important  factors  among  others  that  need  to  be  further  studied  to  properly  design  an  UWB  radar. 
Propagation  models  of  UWB  signals  should  be  developed  and  their  validity  tested  against  actual  data. 
These  models  can  then  be  used  to  determine  the  parameters  necessary  to  develop  viable  detectors.  These 
parameters  can  be  determined  via  appropriate  models  and  compared  by  using  actual  data  to  assess  their 
usefulness. 

The  finding  in  both  Ref.  3  and  in  this  paper  that  UWB  radars  yield  a  processing  gain  over 
conventional  radars  for  targets  which  have  numerous  scatterers  may  imply  that  UWB  radars  have  unique 
advantages. 


7- 
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It  is  recommended  that  future  emphasis  be  placed  on  the  phenomenological  aspects  of  UWB  radar 
signals.  Models  should  be  developed  which  accurately  reflect  the  propagation  of  signals  and  the  effects 
of  clutter  so  that  detectors  can  be  properly  designed. 
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